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Abstract 

Experiments show several steps in hysteresis loops of high spin nanocluster Eeg ■ It is thought that 
these steps are due to thermally assisted resonant tunneling between different quanta spin states. 
Up to now, in calculatiing it, only dipole excitations were considered. Because of the symmetry 
and the power of spin operators in Hamiltonian, we think that other multipole excitation must be 
considered too. In this paper we consider both dipole and quadrupole excitations in Hamiltonian 
and then quantum resonant tunneling phenomena are obtained numerically. As we can show in 
these figures, this phenomenon is different in these two conditions and the second condition is 
nearly to the fact. 
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I. INTRODUCTION 


Today more than five thousand different magnetic molecular clusters are known, and 
about hundred of these show some behaviors that are significant from physics science point 
of view [1]. These molecules are between simple paramagnetic salts and super-paramagnetic 
particles of sub-micrometer size. Thus, they show phenomena such as hysteresis at molecular 
level. Some issues in this subject are similar to those that result in the study of magne¬ 
tization reversal of small magnetic particles, and it is hoped that the molecular systems 
will offer insights into the future, which has applications in magnetic recording and storage 
technologies. 

The molecule Feg (proper chemical formula: [Fe802{0H)i2(tacn)Q\^~^ is magnetic, and 
forms good single crystals.(see Fig.l). In its lowest state, it has a total spin of 10, arising 
from competing antiferromagnetic interactions between the eight S' = 5/2 Fe ions within a 
molecule. Spin-orbit and spin-spin interactions destroy complete rotational invariance, and 
give rise to anisotropy with respect to the crystal lattice directions. A variety of experimen¬ 
tal techniques (electron spin resonance, ac susceptibility, magnetic relaxation, Mossbauer 
spectroscopy, neutron scattering) indicates that their data can be fitted by the following 
model Hamiltonian [2 — 5]: 


n = -DSl - ESl - q^ibS.H (1) 

Where Sx and Sz are the two components of spin operator, D and E are the anisotropy 
constants and are known through a variety of experimental evidences withU = 0.33/c, 
E = 0.092A; [6-7]. The last term of the Hamiltonian describes the Zeman energy which 
is associated with an applied field H. This Hamiltonian defines hard, and easy axes of mag¬ 
netization in X and z direction respectively. 

Magnetization hysteresis are exhibited for this nanocluster below blocking temperature . 
Fig.2 shows typical hysteresis loops with the field applied along the easy axis of magnetiza¬ 
tion. The steps in the loops are clearly visible. These loops become temperature independent 
below 0.35 k and demonstrate quantum tunneling at the lowest levels of energy. 

It suggested that the observed steps in hysteresis loops are due to resonance tunneling 
between quantum spin states. In calculating only dipole excitations were considered [8] and 
we think that other the excitations, such as multipole excitations are related, and must be 
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FIG. 1. Chemical structure in 
nanocluster Fe%. 



- 2mTmin'' 

FIG. 2. hysteresis loops with the 
field applied along the easy axis of 
magnetization. 


considered. 

In this paper, firstly, periodically dependence of spin tnnneling for Fcg in SU(2) gronp 
is calcnlated. In other word, only dipole excitation in Hamiltonian is considered. Dne to 
the symmetry and the power of spin operators in Hamiltonian and being more accnrate, 
other mnltipole excitations mnst be considered. Then, in this stndy we consider both dipole 
and qnadrnpole excitations. Instanton method is nsed to calcnlate tnnneling phenomenon 
between the gronnd states. This method will be described in the next section. 


II. THEORY 

Instantons are classical paths that rnn between degenerating the classical minima of the 
energy. By classical we mean that the path obeys the principle of the least action and 
satishes energy conservation. However, a path along which energy is conserved cannot rnn 
between two minima and can have real coordinates and momenta. Hence, one mnst enlarge 
the notion of a classical path and allow the coordinates and/or momenta to become complex. 

The instanton method is an efficient way of calcnlating tnnnel splitting, both for particles 
[9] and for spin [10]. It is based on evalnating the path integral for a certain propagator 
in the steepest-descent approximation and is designed to be asymptotically correct in the 
semiclassical limit J —)■ cx) or h —0. 

The tnnneling amplitnde in this method is given by [11]: 
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( 2 ) 


A = J2Dke-^" 

k 

The pre-factor Dk results from integrating the Gaussian fluctuations around the k-th 
instanton. For nanocluster Feg, the tunneling amplitude is obtained from the following 
relation [5]: 


A ^ 2Die--^cos{JQ) (3) 

Where A and 0 related to real and imaginary parts of classical action. In this paper, these 
functions were calculated numerically. Because this pre-factor don’t effect on the location 
of quenching points or steps in hysteresis loops, for simplicity we assume 2Di = 1. 

III. CALCULATIONS IN SU(2) GROUP 

We used expectation value of operators Si in SU(2) group to obtain classical energy of 
the system, in other word, in the hrst case only dipole excitations are considered. Then 


5 *+ = A‘f’sin9 
S~ = e~^'^sin6 

= cose (4) 

By placing these values in Hamiltonian (1.1) obtained: 

Ed = —Dcos^e + Esiv?‘9siv?(f) — Acos9 (5) 

Where A = . The minimum of energy is obtained at: 

9 = 7r/2, 0 = 0, TT (6) 

The minima of energy are: 


A,K7r/2,0) =0 

A,K7r/2,7r) =0 (7) 
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By substituting P = cosO , the instanton must satisfy the constraint: 


DP^ -AP + E{1 - P^)sin^(j) = 0 


(8) 


Solutions for the above equation are: 




-A ± A? A AEsi'n?(j){D + EsitPcj)) 


(9) 


2{D + Esiin?(j)) 

The above solutions, despite conserving energy, must satisfy the boundary conditions 
P(0) = P(7r) = 0 . These solutions give rise to what A. Garg called boundary jump 
instantons [3]. The solution whcih satishes the above conditions is: 


P(0) = 


—A + sjA? A AEsin‘^(j){D + Esin‘S(j)) 


2{D + Esin'^(j)) 

If we consider only dipole excitations, A and 0 obtain in the following forms: 


( 10 ) 


A = ReSd = Re P{(j))d(j) 

Jo 

0= r(i-^(0))# (11) 

Jo 


We calculated these functions numerically and substituted in relation (2.2) and tunneling 


amplitude obtained in the form of Fig. 4. If 
held h, the Fig. 3 would be obtain. 



h 


FIG. 3. cos{JQ) versus magnetic 
field h, in this case only dipole ex¬ 
citations is considered. 


plotted function cos(J0) versus magnetic 



FIG. 4. Amplitude of tunneling 
versus magnetic field h, in this case 
only dipole excitations is consid¬ 
ered. 
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This graphic show that there are three quenching points in spin tunneling amplitude, but 
we know from experimented date that the number of these points or steps in hysteresis loop 
are four and for this reason other multipole excitations must be considered. In the next 
section, we will add quadrupole excitation to calculations and obtain these functions and 
graphics. 

IV. CALCULATIONS IN SU(3) GROUP 

If quadrupole excitation is considered, we use spin expectation values in SU(3) group; 


g+ ^ e^'^cos2gsin6 
S~ = e~^‘^cos2gsin9 

= cos2gcos6 (12) 

If placing these quantities in Hamiltonian (1.1), the classical energy would obtain in the 
following form: 


Ed = —Dcos^2gcos^9 + Ecos^2gsin^9sin‘^(l) — Acos2gcos9 (13) 

The minimum of energy obtained at: 


6* = 7r/2,5f = 7r/4, 0 = O,7r 


Also minima of energy are: 


Ecz(7r/2,7r/4,0) = 0 
Ecz(7r/2,57r/4,7r) = 0 


By substituting P = cos9 , the instanton must satisfy the constraint: 


— Dcos^2gP^ + Ecos^2g{l — P‘^)sin^(j) — Acos2gP = 0 
Solutions for the above equation are: 


p{4>,9) 


—A ± + AEcos‘^2gsin‘^(j){D + Esin'^cj)) 

2cos2g{D + Esin'^c/)) 


(14) 


(15) 


(16) 


(17) 
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Similar to the previous section, solution that has the following form is considered: 



FIG. 5. cos(J0) versus magnetic FIG. 6 . Amplitude of tunneling 

field h, in this case both dipole and versus magnetic field h,in this case 

quadrupole excitations are consid- both dipole and quadrupole excita- 
ered. tions are considered. 

If we consider both dipole and quadrupole excitations, A and 0 will obtain in the following 
forms: 


^TT ^7r/4 

A = ReSci = Re / P{4), g)d4)dg 

Jo Jo 

PIT PIT 

0 = / {cos2g — P{(f)^ g))d(f) — / {1 — cos2g)d'y (19) 

Jo Jo 

Where two variables 6 and 0 are related to dipole excitations and g and 7 are related to 
quadrupole excitations. 

We calculated this function numerically and substituted them in relation (2.2). If we 
plotted tunnel amplitude and function cos(J0) versus magnetic field h, quenching points in 
splitting (Fig. 5) or steps in hysteresis loop (Fig. 6 ) would be obtained. 

As we saw in this plot, the number of quenching points for this Nanocluster are equal to 
the number of steps in hysteresis loop. 


V. DISCUSSION 

In this paper, according to the power of spin operators in Hamiltonian, to get the correct 
number of steps in hysterics loop, we used dipole and quadrupole excitations. 
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At the beginning, only dipole excitation was considered, as we saw in Fig. 4, the numbers 
of steps in hysterics loop or quenching points in tunnel splitting oscillations are three and 
this number is less than real one. If we add quadrupole excitation, this number changes to 
four. Fig. 6, and this is equal to the fact. 

Note that, it is possible that the position of these steps or quenching points is different 
from the real state. To be more accurate, we must add another terms as perturbed terms 
to Hamiltonian. 
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